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Influence of surrounding neutral medium on the stability of 


a fluid cylinder carrying an axial current 
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ABSTRACT 


The influence of a surrounding neutral medium on the stability of a compressible fluid 
cylinder carrying a uniform axial current is investigated. It is shown that the presence of 
the neutral medium results in increased effective density of the fluid and thereby leads to 
reduced growth rate for the hydromagnetic instabilities for the fluid cylinder. 


1. Introduction 


The possibility of confining a hot plasma by means of magnetic fields is 
‘extensively investigated [1] both experimentally and theoretically. Since the 
charged particles are “‘tied’’ to the magnetic field lines, it should in principle 
be possible to confine them to a region by a suitable choice of the magnetic 
field configuration. However, the drift motions of the charged particles due to 
the electric and magnetic fields, the magnetic field gradients and the collective 
phenomena like the hydromagnetic instabilities set an upper limit on the time 
in which a plasma can be confined to a given region. 

The essential feature of a magnetically confined plasma is that it is separated 
from the material walls of the container by a region of very low gas pressure. 

Recently Alfvén and Smars [2] pointed out the advantages in surrounding 
a plasma in the presence of magnetic fields by a neutral gas. In particular they 
showed that if a quantity of heat @ is given to a cylindrical plasma per unit 
length by external sources (e.g. by electromagnetic induction) and Q is greater 
than a certain critical value, Q,, one obtains “burn out” of the plasma whose 
temperature increases to very high values, because the thermal conductivity is 
greatly reduced due to the magnetic field. The presence of the neutral gas 
between the plasma and the walls prevents contamination of the plasma with 
the wall material, a fact which is of considerable advantage in reducing the 
radiation losses due to bremsstrahlung. If the density of the neutral gas is 
higher than that of the plasma it is expected that the growth rate of the 
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hydromagnetic instabilities is considerably reduced from its value when th 


model of the plasma by assuming it is a fluid which is compressible and elec ri. 
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ots is surrounded by vacuum. We shall study in this paper a simp 


cally conducting. It will be shown that there is a significant reduction in the 
growth rate of the hydromagnetic instabilities when a cylindrical plasma car. 
rying an axial current is surrounded by a relatively high density neutral medium 


a 2. Formulation of the problem 


In the experiments reported by Alfvén and Smars [2] and Smars and Johans 
son [3], a glass torus with a mean diameter of 10 cm and a bore diameter of 
about 3 cm is filled with helium gas at pressures ranging from 1 to 400 torr. 
The torus of gas, which is made electrically conducting by means of a pre 
ionisation current, forms the secondary winding of a two-turn transformer through 


which a condenser bank of 2.2 uF capacity charged to a voltage of 30 kV is 


discharged. When the discharge is induced inside the high pressure gas, a rela- 
tively small region carries the whole current and this region is surrounded by 
a high density neutral medium. Our aim is to study the stability of the con- 
ducting gas when it carries the electric current and is surrounded by the neutral 
medium. A rigorous theory of the hydromagnetic stability of a toroidal con- 
figuration presents enormous difficulties and therefore we shall only study the 
case of a cylindrical configuration which is simpler to treat analytically. 

To be more specific, we shall imagine an infinitely long cylinder of radius R 
of a compressible fluid of small electrical conductivity carrying a uniform axial 
current which is assumed to be produced. by sources outside the fluid. Let the 
electrical conductivity be so small that the magnetic field is not frozen into the 
medium. (The reason for this assumption will be apparent later.) The conducting 
fluid is assumed to be surrounded by another compressible fluid which is an 
ideal insulator so that no current can leak into it. Let the non-conducting fluid 
extend up to radial distance AR where A>1. Let a cylindrical rigid coaxial 
non-conducting wall bound the non-conducting fluid. 

The configuration will be in equilibrium with the following distributions of 
je the current density, B, the magnetic induction and pe the pressure. (When- 
ever we use the word pressure we shall always understand hereafter the gas 
pressure only.) 


vv (0,9, Jose fon om= Ze. a4 
aloe for R<r<AR, a 
rt daa ry for 0<r<R, (2.2) 
P (0, 3 Uojo R*/r, 0) for R<r< as 

Py=t Uy jo(R?—7*?) for O< <R, 
».-| 2= £ MoJo ( ) . (2.3) 
Po for R<r<AR. 


In eqs. (2.1), (2.2) and (2.3) the suffix e of je, B, and p, denotes that these 
quantities refer to the equilibrium situation; Pp, denotes in eq. (2.3) the pressure 
in the non-conducting medium. Since the fluid is assumed to be compressible 
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the density is not uniform if the temperature is uniform in the conducting 
fluid in which the electric current is flowing. If we assume that the conducting 
and non-conducting fluids have uniform temperatures T, and T, respectively 
and that the equation of state in the simple form 


Ya 


p=nxT, (2.4) 
where » is the total particle density, x is the Boltzmann constant and 7’ is 


the temperature holds, then the distribution of the density, 02, in the equi- 
librium is given by 


o" = (DP, +tojo(R?—r*)), for O0<r<R, 
20 T'y 
any (2.5) 
= ——_ < Ss 
Q2 =~ T. Poy for R<r<AR, 


where m; is the mass of the ion in the conducting medium and m is the mass 
of the neutral atom in the surrounding non-conducting medium. In eq. (2.5) it 
is assumed that the conducting gas is composed of only electrons and singly 
ionised atoms and that the medium is electrically neutral on the average. If 
‘we assume that the nuclei of the ions in the conducting region and of the 
atoms in the non-conducting region are the same, that is, if the plasma is sur- 
rounded by the atomic gas itself, we have m,;~m. Also, in eq. (2.5) the terms 
‘proportional to the ratio of the electron mass to the ion mass are neglected 
and it is assumed that the electron and the ion temperatures are the same. 


3. Basic equations of the problem and the approximations 


If we ignore the transport phenomena like the thermal conduction and the 
viscous dissipation the basic equations of the problem are, in rationlised M.K.S. 
units, 


dv : 
oy iacaalande Maal bone (3.1) 
“2. -oy-y, (3.2) 

i d 

qj (PO) =O (3.3) 
VxXB=poi, (3.4) 
Vv -B=0, (3.5) 
and I =constant. (3.6) 
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where o is the electrical conductivity of the fluid has to be added to the above 
set of equations. However, in the present problem we assume that the electrica : 
conductivity of the fluid is so small that the time scale, tz, of the diffusion of 
the magnetic field lines is much smaller than any other time scale of the prob- 
lem. The relevant time scale with which the we can compare tg is the time 
scale of the hydrodynamic instabilities, which we may call t,. Our assumption 


implies that 
ta/tr< ht (3.8) 


It can be seen from the considerations of order of magnitude of various terms 
in eq. (3.7) that, in view of the inequality (3.8), the second term on the right 
hand side of eq. (3.7) is the dominant term and that the remaining terms can 
be neglected in comparison with it. We can therefore assume that 


v?B=0, (3.9) 


where the right-hand side can differ from zero by terms of the order of the 
ratio tq/Tt,. Using eqs. (3.4), (3.5) and (3.9) we find that 


Vv xj=0. (3.10) 


Eq. (3.10), together with the constraint expressed by eq. (3.6) and the boundary 
condition that the current density vector has no component normal to the 
interface between the conducting and the non-conducting fluids, enables one to de- 
termine the current density distribution inside the conducting fluid uniquely. 

Before preceeding ahead we remark here that the growth rate of the hydro- 
magnetic instabilities for an incompressible fluid cylinder carrying an axial cur- 
rent obtained by Murty [4, 5] on the basis of the approximations of the type 
mentioned above do not differ very much from those obtained by Tayler [6] 
on the assumption of infinitely large electrical conductivity for the fluid, pro- 
vided the magnetic field is produced only by the uniformly distributed electric 
current flowing inside the fluid itself. In view of this fact, we shall continue 
to work with the above-mentioned approximation also in the case of an ionised 
gas which is assumed to be a compressible fluid whose electrical conductivity 
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an arbitrarily small perturbation into a shape given in circular cylindrical 
oording tes with z-axis coinciding with the axis of the fluid cylinder by 


r=R-+a cos @, , . (3.11) 

p= —mb+ kez, (3.12) 
ere m is an integer and & is a positive real number, and a depends on time aoe 
ch that : 
—_ a= exp (wt). (3.13) . 
ee . 


I eq. (3.12) k=2z/A is the wave number, 4 being the wavelength of the per- 
turbation in the z-direction. The amplitude of the perturbation is assumed to 
be small to start with so that — 


a/R <1. 


7A ws 


Tn calculating any physical quantity we shall ignore terms that are proportional 
to the second and the higher powers of a. Further we shall only consider the 
eases when w is a real number and ignore the cases when w is a complex 
quantity. Our aim is to find the dispersion relation between w and k. For any 
given m and &k and other physical quantities if w is found to have positive 
teal value, we conclude that the interface between the conducting and the non- 
conducting fluids is unstable and therefore the equilibrium we have started with 
is unstable for arbitrarily small perturbations of the fluid interface. We shall 
not consider those unstable perturbations that may exist and keep the fluid 
interface unperturbed since these are physically less important in the stability 
considerations than those that do disturb the fluid interface. 

When the fluid interface is disturbed the time-dependent part of all the 
physical quantities is in the first approximation proportional to the quantity a. 
The distribution of the current density and the magnetic field has been dis- 
cussed by Murty [5] in connection with a similar stability problem where the 
fluid cylinder was surrounded by vacuum. Since the electromagnetic quantities 
are determined, entirely by Maxwell’s equations (eqs. (3.4) and (3.5)) and the 
boundary conditions at the boundary of the conducting and the non-conducting 
fluids and independently of the fluid velocities, the presence of a non-conducting 
medium surrounding the conducting medium is immaterial as far as their dis- 
tributions after the perturbation are concerned. In the following paragraph we 
can assume that the perturbed distributions of j and B in eq. (3.1) are known 
and are the same as those obtained in [5]. 

If we take the divergence on both sides of eq. (3.1) and note that the fluid 
was in static equilibrium before the perturbation and. use eqs. (3.4) and (3.10) 
we obtain, neglecting higher order terms, 
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now becomes, after some rearrangement of the terms, 


io ('Teep\! Sat () eee )}=- vient 
AR B) oy ( C +V esr Vp Mod 


where we have further neglected one higher order term and 


+ ~ % \ 
c-(2) oo 


is the acoustic speed. In the conducting medium the acoustic speed, Ox, is 
given by , 


=a)" 


3.18 a) | 
a ( ) 


Ca (>. 


~ 


In the linear approximation we have only to use the zero order (equilibrium) 
distributions of the pressure and the density as given by eqs. (2.3) and (2.5) 
in order to evaluate the second term on the left hand side of eq. (3.17). Using 
eqs. (2.3) and (2.5) we find that in the conducting fluid 


V : 1 ™. a 
Bora-lineleste am 
where r is the unit vector in the radial direction. We see from eq. (3.19) tat 
the terms involving the gradients of the pressure and the density on the left 
hand side of eq. (3.17) drop out if C3, =2%T7,/m; which, however, is not true 
because C;; is given by eq. (3.18 a). The second term on the left hand side of 
eq. (3.17) can be shown, under typical conditions of the experiments at high 


pressure gas discharges, to be smaller than the second term on the right hand 
side of the same equation. In order to do so let us consider the ratio, a, 
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where L is the characteristic linear dimension of the problem and |q| is the 
characteristic growth rate of the instabilities. For a calculation of the order of 
magnitude of the above ratio, the values of the various quantities from the 
experiments mentioned by Alfvén and Smars [2] and Smars and Johansson [3] 
ean be taken as representative. In these experiments L~4x10~? m; m,~9x10-7" 
kg (helium); 7,~+10° °K; the current density ~3x108 Am~*; the mass density 
#4-10-°kgm®* (corresponding to a pre-discharge pressure of 400 torr); the 
value of |w| corresponding to the above mentioned values of the current density 
and the mass density is, according to the results of the present theory, ap- 
proximately 10°s *. When these values are used the ratio, x, turns out to be 
approximately equal to 10~*. We can therefore simplify the eq. (3.17) as 


Rae eal ate La (3.21) 
02, oF Ho? > : 
for the conducting medium. In the non-conducting medium the pressure, p*, is 


given by 


frat tee ie 8) Hh tg (3.21 a) 
Cor ; 


on 


which follows from eq. (3.17) without any further approximations as the current 
density is zero in this medium. The acoustic speed in the non-conducting me- 
dium is given by 


T..\% 
Ca=(r oe : (3.22) 


The superscripts 7 and e for p in eqs. (3.21) and (3.21 a) are used to distinguish 
the pressure in the conducting and the non-conducting media respectively. 
4. Pressure distribution 


In eq. (3.21) the current density distribution is given [5] up to first order 
terms by 


i j= (rs Jo, Jz)> (4.1 a) 
h jp = —ak Tn (lt) sin 
where it jo Tr ER) 
. Sami, Iglkr) (4.1) 
ese TER) 
atthe Lather) 
and I2=Jo akjyr (é R) COS @, 


> 


fap “Sd Bi 6) Online halen Se ROBE 
Vv? ph (r) = — Moi, rer aden, ise J 
i BALSA! Re IE Phy Ta Oa ae a 


5 
7 bast 
c ee’ 


ai 
ae 


MiLETR 


re Maen pine 
and (- Ze Sat Vv api (r, 0, z) * Pl (eR) COS Y. , 
: =_ ing ‘3; 21 : 
Since a was assumed to be given by eq. (3.13), eq. (4.5) becomes i a 
Ue 
2 o In (Er) ae 
—2_ +9?) ol (r, 0:2) =2'u, ke 2 =~ cos 9, (4.6 
( oat V ) rier » 2) ROT a MO? 
which suggests that pji(r, 0, z) can be separated in the variables as : dha 
pi(r, 0, z)=pi(r) cos g, = (4.7 
so that pi (r) satisfies the equation 
led d 5 me ie, Gi of flex) 4.8) 
EF a (rar) — (7a) | ho 20 0k hy ae 
w. 
where P= (3. + #) ; (4.9) 
Cs1 


It can be shown that the general solution of the inhomogeneous ordinary dif- 
ferential eq. (4.8) is given by | 


2 yk 76 In (ker) 
Epa 0 
pi (r) =A, I, (br) + e—P Ti, (k RY’ (4.10) 


where A, is an arbitrary constant and J m(lr) is the modified Bessel function 
(regular at the origin) in whose argument the quantity J is defined as the posi-_ 
tive square root of the right-hand side of eq. (4.9). In the non-conducting fluid 
the radial part of the first order pressure distribution is given by ; 
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; of Re fluid’ eters is zero Biers: Since iat are no elect. 

orces in this medium, this condition is equivalent to the vanishin; 

the na derivative of the first order pressure r= AR, which means (see eq. 
11)) that 


_Lin(AhR) + Dy Kin (Ab R)=0, 


where the prime to the Bessel functions denotes differentiation with respect to 
the sen. The continuity of the pressure at the deformed boundary avo. 


by eq. (3.11) gives another condition. The pressure in the conducting fluid, p', 
is given by eqs. (4.3), (4.7), (4.10) and (4.13): ’ 
9 
pi=E—kuyjor* taf, In (UR) + 5a ae Jn ie COs @. (5.2) 


=T he pressure in the ee ae fluid is given by 

g "pe =F +a {D, In (hr) + D, Kn (hr)} cos @, (5.3) 
which is the general solution of eq. (3.21a), F being an arbitrary constant. The 
continuity of the pressure at r=R-+a cos p means that p' equal p° there. 
-Equating the terms of the corresponding order we obtain 


Z 


E-}ujR =F, (5.4) 
3 92 Im (R ; 
and yids eee ek Tn (UR) =Dyin(hR)+D,Kmn(hR), (5.5) 


where we have neglected higher order terms proportional to a’. In eq. (5.4) F 
can be identified with p, in eqs. (2.5) and therefore E is determined. Here- 
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after we need consider only eq. (5.5) in connection with obtaining the dispersior 
relation. 

In addition to the two eqs. (5.1) and (5.5) we can obtain as shown in the 
following two more equations for the four unknowns oA; D, andeD;. 
virtue of eqs. (3.11) and (3.13), eq. (3.1) can be rewritten as 


j | 

054" cos p-{-22+GxB)] . (6.6), 

6 , 

: 

where the suffix ‘b’ to the brace and.in g, indicates that the respective quan- 

tities must be evaluated at boundary approached from either side. The only 

quantity in eq. (5.6) that is not listed so far is the distribution of B. As was 

mentioned in §3 we can write down at once the distribution of B in the present 

problem from [5] where it is calculated for a case when there is an infinitely 

conducting coaxial sheath at some distance from the conducting fluid. For the 

present calculations we need only the distribution of B in the absence of such 
a conducting sheath and it is given by (note the change in units) 


B=(B,, Boe, B.), (5.7 a) 
where B,= {4 In (kr) eee eas sin g, 
nimi fat ? rnaS n| 
and B,=AI,(kr) cos 9, 
where A= —pUyamiy Jn (kB) Bn (kB) (5.7 b) 


I, (k R) : 


in the region 0<r<R-+a cos y. We shall not write the distribution of B in the 
outer region as it is not required for our purpose. If we substitute for o LP 
and B from eqs. (2.5), (4.1), (5.2) and (5.7) respectively in eq. (5.6), we obtain, 
using eq. (5.7 b) 


2 , ; 2?) KRIn(kKR) m?1,(kR) 
w°=—-A,lI,(LR)+14 slitz im —E m 
@ in (2) Hoi | Fae to 7 eR) ERI. (ER) 1 R) Kn (ER), 
(5.8) 
mM; 
where Ge) Po (5.9) 


Similarly if we substitute for p® from eq. (5.3) and 0=@., the density of the 
outer medium in eq. (5.6) we obtain 


0,0" = — D, hI, (hR)— Dh K;, (A R). (5.10) 


Eliminating the constants A,, D, and D, among th f 
and (5.10) we obtain il : g the tour eqs. (5.1), (5.5), (5.8) 
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oO for given values of the remaining parameters. 


eee oS 6. Discussion 


We shall discuss eq. (5.11) in the following limiting cases, St 


¥ 


(a) Both fluids incompressible 


Z It we assume that both the conducting and the non-conducting fluids are in- 
vpressible, that is, Csi and C.2 have infinitely large values, we can show 
t eq. (5.11) reduces to : 


~ 


wo" [0, +02 Q(m, A, KR)]=poj0P (m, kR), (6.1) 


Tin (RB) Im (BR) Kin (A&R) —In(AKR) K(k R) 


: = = Sa PT Te 6.2 
where — Q(m, A, kR)=7 ER) In (ALR) Kn (ER) ~ Kn (ABR) In (ER) oe) 
a _ | BR [Ln(kB) In aa 

— Rea Ag: ree) I, (kR) 

E m? In (kR) 

J —_ m(kR) K(k R)S. 6.3 
If we further assume that the conducting fluid is surrounded by vacuum, that 
is, if 9,=0 in eq. (6.1), the resulting equation! is the same as the one derived 


and discussed in detail in [4,5] where it was shown that w?>0 for values of 


ie Sree 


, 
i 


TRA Ce Sa 


mn In this ‘connection see also the appendix at the end of this paper. 
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Fig. 1. Fig. 2. 
Fig. 1. Variation of Q(1, A, x) with x=kR for A=2 and oo. The function Q is a m 
of the effective density of the incompressible and non-conducting fluid surrounding th 
compressible and conducting fluid cylinder for a perturbation of a given wavelength A=: 
f is the radius of the conducting fluid cylinder. The value 1 in the arguments of the fun 
Q indicates that Q corresponds to the helical perturbations of the fluid cylinder. A is the 
of the radius of the rigid and non-conducting coaxial cylinder to R. As x tends to infinit 
large values, the value of Q(1, A, x) tends to 1 for both A=2 and oo, 
Fig. 2. Variation of * (in non-dimensional units) as given by eq. (6.1) for helical pertur! 
tions with x=kR, for two values of A, the ratio of the radius of the rigid and non 
ducting coaxial cylindrical wall to R. The influence of this cy'indrical wall is negli 
particularly for perturbations of the fluid cylinder with very small values of the wavele: 
along the axis. The ratio Q,/0, is assumed to be 10. 


~ 


kR>0 (which is the only range of the values of kR which is of interest t 
us) when |m|=0 or 1. It was also shown there that when |m|>1 there exis 
a range of values of kR from zero up to a certain value, (kR),, such tha’ 
w* 20 according as kR=(kR),. Since w?>O0 corresponds to the unstable mode; 
of the perturbation the fluid cylinder was shown to be unstable for all th 
values of m and the growth rate was shown to tend to the same maximu 
value for all the values of m as the value of kR tends to infinitely large values 
If e.+0, only the left hand side of eq. (6.1) is effected while the right hand 
side is unaltered. The function Q(1, A, kR), which may be considered as 
measure of the effective density of the non-conducting medium for a perturba- 
tion of a given wavelength, is plotted versus &R in Fig. 1 for two values of A, 
namely, 2 and co. It is seen in this figure that the function Q is positive for 
all the positive values of kR and that when kR>3 the value of Q(1, 2, kB 
is nearly the same as that of Q(1, co, kR). Both the functions will approach 
the limiting value 1 as kR tends to infinitely large values. For the same values. 
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ch shows that the growth rate of the hydromagnetic instabilities of an in- 
pressible fluid cylinder carrying an axial current is significantly reduced 
en the fluid cylinder is surrounded by a non-conducting and incompressible 
d of relatively high density. 


ar 


(b) Incompressible conducting fluid surrounded by compressible 
non-conducting fluid 


_ In the high pressure gas discharges where the current carrying channel is sur- 
rounded by a non-conducting gas, the temperature of the conducting gas is larger 
than that of the non-conducting gas by a factor of hundred or even more and 
therefore the conducting gas has a correspondingly higher acoustic speed than 
ee in the non-conducting gas. Under such conditions an approximate value of 


he growth rate of the instability can be obtained by taking the ideal case of 
an incompressible conducting fluid surrounded by a non-conducting compressible 
fluid. The dispersion relation under these conditions can be shown to be given by 
A 


kIi,(kR : 
wo Jere: re aieipdes A, ms) = [ojo P (m, k RB), (6.5) 


which is obtained from eq. (5.11) by letting C1 tend to infinitely large values. 

A complete solution of eq. (6.5) for w as a function of kR for given values 
of the remaining parameters is not simple because the equation is an implicit 
equation for w, but we notice that both sides of the equation involve only real 
functions of real arguments if w is a real quantity which shows that there can 
be real values of w which are solutions of eq. (6.5). Among the values of a, 
those corresponding to large values of & are of particular interest since the in- 
stabilities have in general larger growth rates, the larger the value of & is. An 
examination of eq. (6.5) shows that w” cannot have infinitely large value for 
any value of kR. Therefore we see that 


w"/k? Cs < 1, 
when &£R—-co, that is, w tends to a finite limit when k tends to infinitely 
large values. If we expand all the functions in eq. (6.5) in a power series in 


w?/k? C2, and 1/kR and retain only the leading terms, we find that w? has a 
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finite value plus higher order terms which vanish in the limit of infinitely larg 
values of Y The limiting value of |w| thus obtained will be found to be t 
same as that when both the fluids are assumed to be incompressible. This mean 
that the growth rate of perturbations of infinitesimally small wavelength in the 
present case also is given by eq. (6.4). 


(c) Both the fluids compressible 


If both the conducting and the non-conducting fluids are assumed to be com- 
pressible and the non-conducting fluid is bounded by a rigid coaxial cylindrical 
wall, the dispersion relation is given by eq. (5.11) where both the quantities l 
and h occur. We see that this dispersion relation also involves only functio 
that are real functions of a real argument if m is a real quantity. Followi 
the arguments that are employed in the previous paragraph, one can conclud 
that the limiting value of w” obtained in this case when kR tends to infinitely 
large values is also the same as the one obtained when both the fluids are as- 
sumed to be incompressible and that the growth rate of the perturbations of 
infinitesmally small wavelengths is given by eq. (6.4). 


ve Concluding remarks 


We have shown that a conducting fluid cylinder carrying an axial electric 
current and surrounded by a non-conducting fluid which is bounded at some 
distance from the axis is less unstable than when it is surrounded only by 
vacuum, the remaining conditions being the same in the two situations. The in- 
fluence of the surrounding neutral medium is manifested only in the increased 
effective density of the conducting fluid. The effects due to the compressibility 
of the fluids can be ignored for perturbations of infinitesimally small wavelength. 
The growth rate of the instabilities has a maximum value for perturbations of 
vanishingly small wavelength. 

An important physical factor that should be taken into accunt when the wave- 
length of the perturbation becomes very small is the viscosity of the fluids. 
Even though it is very difficult to estimate the viscosity of the ionised gas in 
the presence of an electric current, its non-zero value, however small it may be, 
will lead to viscous stresses that cannot be ignored for perturbations of very 
small wavelengths. As we have already seen, the perturbations of vanishingly 
small wavelengths have the largest growth rate in the absence of viscous forces. 
Due to the damping effects introduced by the viscous forces one can expect 
that the growth rate of the unstable modes should tend to zero for perturba- 
tions of vanishingly small wavelength [6]. This gives rise to the result that there 
is a wavelength of the perturbation for which the perturbation grows at maxi- 
mum rate in the presence of viscous forces. We can obtain an idea of the ap- 
proximate value of the wavelength of the perturbation that grows faster than 
the perturbations with other values of the wavelength from the calculations of 
Tayler. He showed that a compressible and ideally conducting fluid cylinder 
with surface current distribution and surrounded only by vacuum has a maxi- 
mum growth rate corresponding to a perturbation of a certain wavelength, Am, 
whose approximate value is given by 
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= Am =6 RB (d/R)S, (7.1) 
 . 

where d is the mean free path of the particles in the compressible fluid. The 
value of /,, is in general smaller than R since d < R, a condition which is im- 
plicity assumed to be valid when the fluid ‘approximation is made in stability 
calculations. The value of the wavelength of the perturbation that grows at a 
maximum rate in the present problem also can reasonably be expected, although 
we have not proved it, to be of the same order as Am given above. Therefore, 
the growth rate of the perturbations of infinitesmally small wavelength, given 
by eq. (6.4), can be taken to be a first approximation to the actual growth 
rate when viscous forces are taken into account. We can only say from our 
previous calculations that the growth rate obtained for perturbations of van- 
ishingly small wavelength is an upper limit and the actual growth rate is smaller 
than the one given by eq. (6.4). 
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Note added in the proof: In the references (4) and (5) given in the present paper the name 
“Taylor” should be read, whenever it occurs, as ““Tayler’’. 


APPENDIX 


In [4,5] Murty has obtained a dispersion relation when e,=0 taking into 
account the influence of the surface tension of the fluid. It can be shown that, 
when 0,+0, the.required dispersion relation for the case of incompressible fluids 
including the effect of interfacial tension can be obtained by adding to the right 
hand side of eq. (6.1) of the present paper a term 

ThRI,, (kB) 


a (FR? +m? —1), 
Bi, (kR) | is 


where 7 is the interfacial tension between the conducting and the non-conducting 


fluids. 
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